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The performance of the GlueX Forward Calorimeter was studied using a small version of the detector and a variable energy electron 
beam derived from the Hall B tagger at Jefferson Lab. For electron energies from 110 MeV to 260 MeV, which are near the lower- 
limits of the design sensitivity, the fractional energy resolution was measured to range from 20% to 14%, which meets the design 
Q^goals. The use of custom 250 MHz flash ADCs for readout allowed precise measurements of signal arrival times. The detector 
^*<( achieved timing resolutions of 0.38 ns for a single 100 mV pulse, which will allow timing discrimination of photon beam bunches 
^ and out-of-time background during the operation of the GlueX detector. 
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1. Introduction 

The upcoming GlueX Experiment [1] in Hall D at the 
Thomas Jefferson National Accelerator Facility (JLab) is an 
experiment that will primarily search for mesons with exotic 
quantum numbers. Meson states of interest will be produced 
using a photon beam incident on a proton target. It is crucial 
for the experiment to have a detector with good resolution and 
high acceptance for multi-particle events - studying the angular 
distributions of decay products is an essential tool in identify- 
ing underlying structure of the produced mesons. The Forward 
Calorimeter (FCAL) will be an essential detector for the ex- 
periment, providing energy measurements and timing informa- 
tion for photon showers in the forward region with polar angles 
6 <\2° and energies between 100 MeV and 5 GeV. 

In the fall of 2011, a 25 -element miniature version of the 
FCAL was constructed for a beam test underneath the exist- 
ing photon tagger of Hall B [2] at JLab. The main goals were 
to verify, as expected, that the hardware configuration of the 
FCAL modules would meet the desired energy and timing res- 
olutions. While many studies with prototypes were performed 
during the design phase, this represented the first test of the pro- 
duction hardware and data acquisition system in a beam. Below 
we discuss the results of this beam test and compare them with 
previous measurements. 



2. Setup of Experiment 

2.1. FCAL components 

The FCAL for the GlueX experiment consists of 2,800 lead 
glass modules, each coupled to its own type FEU 84-3 photo- 
multiplier tube (PMT) and Cockcroft-Walton base (similar to 
that detailed in Ref [3]). The lead glass blocks were equivalent 



to type F8 manufactured by the Lytkarino Optical Glass Fac- 
tory [4], and each have transverse dimensions of 4 x 4 cm^ and 
are 45 cm long. The Cherenkov light emitted by the electro- 
magnetic showers produced within the lead glass blocks will be 
detected by the PMTs. The resulting PMT current pulses are 
digitized by 12-bit 250 MHz flash analog-to-digital converters 
(fADCs) designed by JLab [5]. Figure 1 shows an expanded 
view of one of the 2,800 FCAL modules. 
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Figure 1: Expanded view of a single module of the GlueX 
FCAL. 

The lead glass blocks and most of the PMTs are the same 
as those used in previous experiments, E852 at Brookhaven 
National Laboratory [6, 7] and the RadPhi Experiment at 
JLab [8, 9]. 

In the GlueX experiment, magnetic shielding of the PMTs is 
necessary due to the stray field of the solenoid magnet of up 
to 200 G. The PMTs are placed well within a double layer of 
soft iron and mu-metal shield. The light collection from the 
lead glass block to the PMT is facilitated using a cylindrical 
acrylic light guide glued to the PMT. An Eljen EJ-560 optical 
interface "cookie" is used to connect the guide to the lead glass 
block. The primary goal for this beam test is to verify that we 
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can achieve the design energy resolution at low energies, where 
statistical fluctuations in the number of Cherenkov photons de- 
tected dominate the resolution. This, in turn, would validate 
the design of the light collection system using optical photons 
with the characteristic Cherenkov angular distribution inside of 
the block, a task that had only been performed with ray-tracing 
simulation. 

Besides measuring the energy of electromagnetic showers, 
the FCAL will measure the timing of the showers. During 
GlueX running, the JLab Continuous Electron Beam Accel- 
erator Facility (CEBAF) beam will provide an electron beam 
bunch every 2 ns. Electrons in the bunch interact with a thin 
diamond crystal radiator to produce a linearly-polarized photon 
beam. During nominal running configuration the diamond is 
oriented so that the coherent bremsstrahlung process produces 
a peak in both linear polarization and photon flux at photon 
energies of about 9 GeV. The photon energy is measured by 
measuring the momentum of the recoil electron using the Hall 
D tagger. The photon rates are expected to be up to 10*^ s ' in 
the coherent peak between 8.4 and 9.0 GeV, leading to a high 
rate of signal photons as well as electromagnetic background. 
The timing measurements of the FCAL will allow identification 
of the respective beam bunch that created the photon of interest 
and also help in reducing accidental backgrounds. 

2.2. Setup of Detector 

A 5 X 5 array of FCAL modules was constructed at Indi- 
ana University. The modules were encased in a light-tight alu- 
minum body, with a gear drive that allowed the modules to be 
tilted at an angle so that the front face of the modules were per- 
pendicular to the incoming electron trajectories. Figure 2 shows 
a drawing of the detector 
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Figure 2: Drawing of the array from the front and back with 
the trigger box attached. For illustration purposes, the lid of the 
trigger box is shown open. 

For triggering purposes, two arrays of five Eljen EJ-200 plas- 
tic scintillator paddles were placed at the front of the detector. 
The three inside paddles had widths of 4 cm, while the two out- 
side paddles were 6 cm wide and used as veto signals to the 
trigger. The two arrays were oriented perpendicularly to each 



other, with the centers of each paddle centered on the border 
of two modules, as shown in Fig. 3. The scintillators spanned 
the length of the perpendicular array of modules. When look- 
ing into the detector from the upstream end (the "beam's eye 
view"), the horizontal paddles were labeled hl-h5 from the top, 
and the vertical paddles were labeled vl-v5 from the left. Each 
of the trigger paddles were optically connected to the same type 
FEU 84-3 PMTs used for reading out the detector modules. The 
signals from the PMTs were converted into NIM signals by a 
discriminator, and the digital output of the discriminators were 
recorded into the data stream using the same fADCs that were 
used for recording the pulse height of the lead glass modules. 
The trigger arrays were placed in a dark box that was attached to 
the front of the main detector. The dark box was made of 1 .27 
cm-thick plywood and painted black, with a thin (L27 mm) 
aluminum plate underneath the trigger setup to hold the trigger 
paddles in place. 
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Figure 3: The placement of the trigger scintillator paddles on 
top of the modules. The position of the front face of the 5x5 
array of detector modules is shown as the dashed lines. 

In addition to the main trigger counters, there was a thin re- 
mote paddle located several meters upstream of the detector, 
along the nominal trajectory of electrons hitting the detector 
(shown in Fig. 5). The remote paddle was made of a 0.125 
cm-thick plastic scintillator, wrapped in aluminized mylar, and 
connected to a Photonis XP2020 PMT. The signal was also con- 
verted into a NIM pulse signal using a discriminator and read 
out with a fADC channel. The purpose of this remote paddle 
was to identify electrons following the nominal electron trajec- 
tory during off'-line analysis. 

Figure 4 shows a diagram of our trigger logic, which was 
handled by standard NIM modules. The PMT pulses from each 
of the 25 modules within the detector array were digitized using 
a fADC. To trigger the readout of an event event to disk, we 
required a hit within at least one of the inside paddles from 
each array. The outside paddles were used as vetoes, and the 
veto signal was also individually recorded. 

2.5. Test Setup in Hall B 

The beam test was conducted within Hall B at JLab, which 
housed the CLAS detector and its photon tagger setup. Dur- 
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Figure 4: Electronics diagram for the trigger. 



ing our beam test, the CLAS Collaboration was running the 
HD-ICE experiment with a real photon beam produced by 
bremsstrahlung, allowing us to use the electrons that had radi- 
ated photons. Our beam test was conducted parasitically to this 
run, and therefore we had no control over the beam energy or 
intensity. Figure 5 shows a schematic of Hall B and the tagger 
setup. The electron beam from the accelerator comes in from 
the left and is incident on a very thin radiator foil. Immedi- 
ately following the radiator is the CLAS tagger magnet, which 
is a uniform-field dipole magnet. The electron trajectories are 
determined by the energy of the electron, and the position and 
angle at the exit of the magnet can be calculated. Details of the 
CLAS tagger setup can be found in Ref. [10]. 




Figure 5; Schematic of the Hall B and the CLAS tagger sys- 
tem. The electron beam radiates a photon at coordinates (0,0). 
The trajectories of degraded electrons with varying energies are 
given by the lines originating from the vacuum exit window. 
Depictions of our detector are shown in the alcove area and 
floor. 

The electrons exit the vacuum of the tagger magnet through 
a very thin exit window [11] and pass through the CLAS tagger 
detectors called the E-counters and T-counters. The E-counters 



are a plane of 384, 4 mm-thick plastic scintillators that detect 
the position of each electron for energy identification. The T- 
counters form a second plane of 2 cm-thick plastic scintillators 
that are used for timing identification of the electrons. The two 
detector planes are separated by 20 cm, and each counter within 
each plane is rotated so that it is facing normal to the incoming 
trajectory of electrons. 

In our beam test, the detector was placed mainly in the area 
seen in the left side of Fig. 5 that we call the alcove. The detec- 
tor was placed at the calculated location that intercepted elec- 
trons with 5% of the accelerator beam energy and was angled 
to face the incoming electrons. We also placed the detector in 
the floor area (right side of Fig. 5), where data were collected at 
a higher electron energy fraction. Due to mechanical interfer- 
ence and running constraints, we could not properly align the 
detector in the floor area, and hence did not use those runs for 
our energy resolution measurements (Sec. 3). However, since 
the timing resolution studies (Sec. 4) are not very sensitive to 
the detector alignment, we utilized these runs for this purpose, 
as they had the largest range of signal sizes. 

3. Energy Resolution Study 

3.1. Data Structure 

For the energy resolution studies, our detector position was 
fixed and data were taken at three separate accelerator energies, 
allowing us to measure the energy resolution for these energies. 
The three runs that were used for our results had accelerator 
energies of 2.537, 4.446, and 5.542 GeV. Each run had approx- 
imately 2 X 10^ events. 

Each channel of the f ADCs samples the PMT output voltage 
every 4 ns, and 50 samples were stored for each event. Figure 6 
shows the read out samples from a single channel in one event. 
However, zero-suppression was implemented by setting a pre- 
programmed threshold and only channels that had at least one 
signal sample larger than this threshold were written to the out- 
put stream. The threshold applied to each channel was based on 
the pedestal, or baseline, which was measured without any in- 
put signals. The threshold was set to 10 ADC counts above the 
pedestal. Off-hne studies showed that the resolution worsens 
with larger thresholds. After data taking, it was discovered that 
the thresholds for modules in the top two rows were inadver- 
tently set much higher, leading to a significantly worse energy 
resolution. Therefore, the sub-sample of events with showers 
centered on these modules were not used in the determination 
of the energy resolution (see Sec. 3.4). 

The module response to energy deposited in the lead glass 
is obtained by summing the pedestal-subtracted samples over a 
range of time. This range was optimized and was 20 samples 
(or 80 ns) starting at sample 9. The sum was performed ofi"- 
line for the beam test, although for the GlueX experiment the 
determination of the signal size and time will be computed in 
Field-Programmable Gate Arrays (FPGAs) to manage the size 
of the output data stream. 

The first step in our analysis was to ensure that we had a sam- 
ple of events that had one and only one electron signal with the 
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Figure 6: ADC values read out for one channel in one event. 
The pedestal which is given by the average of the first eight 
values is shown as the dashed red line. A clear signal is seen 
for this event. 



correct timing. Therefore, a preliminary skim was applied to 
the data to extract such events. For each inside trigger scintilla- 
tor, we determined the pedestal by taking the average of the first 
eight samples recorded. To exclude events that had an unstable 
pedestal, we removed events that had channels with a pedestal 
RMS greater than 2. This eliminated events that were on the tail 
of a previous event and also noise fluctuations at the beginning 
of the event. 

Next, to ensure that the timing of the event was correct, the 
recorded NIM signals were examined and a coincidence in the 
leading edges was required. Also, events with a prolonged or 
second signal were removed. Finally, we required that only one 
horizontal and vertical trigger combination satisfied all of our 
trigger conditions to ensure that there were not two separate 
incoming electrons. 

Similar requirements on the pedestal RMS and timing of the 
signal were made on the remote paddle based on the distribu- 
tion of the data. We verified no signal was present on the veto 
just before or after the triggered events. For all detector mod- 
ules, if any sample indicated an underflow or overflow in the 
ADC counts, that event was removed. The requirements on the 
trigger counters removed 10-20% of the data, while the require- 
ments on the remote paddle removed 70-80% of our initial data 
sample indicating that our detector was frequently triggered by 
electrons that didn't pass through narrow angular range that was 
spanned by the remote paddle at the exit of the tagger window. 

3.2. Gain Balancing Procedure 

The total signal size for the event is given by summing the 
signals from all modules. For each run, the distribution of 
the total signal is fit with a Gaussian function and linear back- 
ground. The width (cr) of this distribution divided by the cen- 
troid is then the measured resolution. Below we outline the 
procedure that was used to optimize the resolution consistently 
over all runs. 

Prior to the beam test, the gain characteristics for the PMTs 
for each module were measured, and the high voltage (HV) set- 
tings were adjusted to equalize the gains. With these initial 
settings, data were taken to measure the initial energy resolu- 
tion. A software gain balancing technique was then applied to 



minimize the event-to-event variance of the sum over all mod- 
ules by scaling the signal size from each module by a constant 
module-dependent gain factor This minimization can be for- 
mulated in terms of a Lagrangian multiplier method, which is 
described in detail in Ref. [8]. In Ref. [8] the mass squared of 
the tt" and rj mesons were used as constraints and measured dif- 
ferences from the meson masses squared were minimized. We 
fixed the average total signal size (in arbitrary units), and the 
standard deviation from this average was minimized. Such a 
problem is linear with respect to the gain factors and can be 
solved exactly, so that no iteration of the minimization is nec- 
essary. Next, using the computed gain factors, we adjusted the 
HV values based on knowledge of the gain characteristics of 
each PMT to better equalize the gains of the PMTs. Finally, we 
repeated the gain factor determination, and utilized these gain 
factors in our analysis. 

We explored the stability of our software gain-balancing pro- 
cedure by selecting events whose sum were within +3cr or ±5cr 
of the mean and determined gain constants for both sets. This 
selection reduces the influence of non-Gaussian tails on the 
computation of the variance and the gain constants. Our re- 
sults are therefore given as the average resolution from these 
selections, with an error given by half the diff'erence, which is 
negligible on the scale of other systematic uncertainties in the 
determination of the energy resolution. 

Figure 7 demonstrates the method of determining the reso- 
lution and shows the distribution of the total signal from all 
modules for two given runs at the same incoming electron en- 
ergy. Figure 7(a) shows a run in which we used our initial HV 
settings without the PMT HV adjustments, and Fig. 7(b) is for 
a run when it was applied. For each case, the black hollow 
histogram is the distribution before the software gain balanc- 
ing procedure is applied, and the gray (red online) hatched his- 
togram is after We see that software gain balancing is essential 
to obtaining optimal resolution, and the resolution is best when 
software corrections are applied to approximately equal-gain 
modules to start, i.e., the precision that is lost when a low gain 
module is digitized by the fADC cannot be recovered by a sim- 
ple multiplicative constant. In principle, one can continue to 
iterate hardware HV adjustments based on software gain con- 
stants; however, we performed only one such HV adjustment 
and utilized the software gain balancing procedure for all sub- 
sequently analyzed data. 

3.3. GEANT Simulations 

Corrections for energy loss and scattering in the material be- 
fore the electrons entered the detector were determined using a 
GEANT4-based [12] simulation program. For each simulated 
data set, nominal energies for the electrons were chosen to cor- 
respond to actual electron energies used in the data collection. 
The fractional spread of the actual CEBAF electron beam is of 
order 0.01% and is negligible. Electrons were produced with 
energies centered around this value and with a range of ±10% 
to cover the range of rays that would scatter into our acceptance. 
For each event, the incident angle and position with which the 
electron exited the vacuum window were calculated based on 
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Figure 7: The distribution of the total signal of all modules be- 
fore (hollow/black online) and after (hatched/red online) soft- 
ware gain balancing. Fits with a Gaussian function and Unear 
background are shown by dashed lines. The two plots are from 
runs with the same electron energy; the top (bottom) plot is 
before (after) PMT HV adjustments. For the bottom plot, in 
addition to balancing gains, the overall average HV settings for 
the PMTs were raised, which increased the average total signal. 



the CLAS tagger geometry. The electron energy distribution 
was modeled on a realistic bremsstrahlung distribution. 

Using the simulation, we recorded the distribution of total in- 
coming energy into our detector. The same trigger conditions 
as used in the analysis of the data were imposed on the sim- 
ulated events. We fit the incoming energy distribution for our 
final selection of events with a Gaussian function to determine 
the mean and width. We concluded that the mean incoming en- 
ergy into our detector was 12-14 MeV lower than the nominal 
energy with a spread of 4-6 MeV. Table 1 summarizes the re- 
sults of our simulation for all energies. With the assumption 
that the intrinsic detector resolutions and the spread of energies 
add in quadrature, we can subtract this contribution to obtain 
the intrinsic detector resolution. 

Table 1: Predicted energy and spread of electrons incident on 
the detector during three run periods. Eq is the nominal incom- 
ing energy, while E and crioss are the center and width of the 
distribution obtained from the simulation, respectively. 

Eo (MeV) E (MeV) en,,, (MeV) crioJE (%) 



126.85 


114.19 + 0.03 


4.54 ± 0.02 


3.98 


222.30 


209.29 + 0.05 


5.15 + 0.03 


2.46 


277.10 


263.70 ± 0.04 


5.70 ± 0.03 


2.16 



is the quadratic sum of the intrinsic detector resolution, crjet, 
and the spread due to different incoming energies, which we 
approximate with our simulation results, crioss/E. With the as- 
sumption that the intrinsic detector resolution and the spread of 
energies adds in quadrature, we can remove the spread of ener- 
gies to obtain the intrinsic detector resolution. 

3.4. Energy Resolution 

Figure 8 shows the energy resolution for five different trig- 
ger combinations. For all runs, the statistical uncertainties from 
the fits to determine o-jj,f.^s were negligibly small compared to 
the systematic variations observed by changing the analysis 
method. The dominant source of systematic uncertainty is vari- 
ation in the performance of the individual modules. The nine 
different trigger possibilities, 3 horizontal x 3 vertical, could be 
analyzed independently and each trigger combination populates 
the individual modules differently. Therefore, the energy reso- 
lution varies with the location of the incident electron. Three 
trigger combinations were excluded due to the incorrect hard- 
ware threshold settings and a fourth combination was also ex- 
cluded because it consistently gave poor resolution, although a 
clear reason was not identified. Each marker represents a dif- 
ferent trigger combination, and have been shifted horizontally 
slightly for visual clarity. The measured energy resolution for 
the three electron energies of 1 14.2+4.5 MeV, 209.3+5.2 MeV, 
and 263.7 + 5.7 MeV, as given by the average and standard 
deviation of the five trigger combinations, are 19.8 ± 0.5%, 
14.3 ± 0.4%, and 14.1 + 0.5%, respectively. 

We expect that the dominant contribution to the resolution 
in this energy regime is statistical fluctuations in the number of 
photoelectrons produced in each PMT for each electron shower. 
Based on simulations of the light collection optics that were 
carried out in the design phase, the expected resolution of the 
GlueX FCAL can be parametrized as 
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We incorporate these simulation results into our analysis in 
the following way. The spread in measured energies, o-jjif^-^s. 



This is shown as the solid blue curve in Fig. 8. Our energy res- 
olution results are consistent with the design goal in this energy 
regime. 

4. Timing Resolution Study 

The second important performance metric of the calorimeter 
is the resolution of the absolute anival time of the PMT pulse 
that is digitized by the fADCs. The fADCs for the FCAL sam- 
ple the PMT pulses only every 4 ns; however, it is possible to 
utilize several samples and knowledge about the pulse shape to 
determine a pulse time with resolution much better than 4 ns. 
In a previous article [13], the timing resolution of the fADCs 
was tested using a pulsed LED source, and a timing resolution 
of less than 1 ns was determined for pulse heights of 100 mV 
and higher. In the following, we apply the same methods as in 
Ref. [13] for the purpose of confirming that we can obtain ad- 
equate timing resolution using actual electromagnetic showers. 
For this analysis, we utilize a single run with 507,828 events 




0.10 



0.15 0.20 0.25 

energy at detector (GeV) 



Figure 8: The energy resolution results against incoming en- 
ergy. The different markers represent different trigger combi- 
nation selections, and are slightly horizontally shifted for vi- 
sual clarity. The solid blue curve is the expected resolution for 
GlueX. 



that was taken at an incoming electron energy of approximately 
1275 MeV. The timing resolution depends strongly on the PMT 
signal amplitude, and the runs at this energy allowed us to mea- 
sure the timing resolution across the largest range of sample 
pulse heights. 

4.1. Timing resolution method 

An example of the digitized PMT signal by the fADCs was 
shown in Fig. 6. We use the time at which the pulse reaches 
half the recorded maximum to specify the signal arrival time. 
This time, which we call t^, is chosen because it is in the region 
of maximum slope. 

The timing resolution was determined using a linear inter- 
polation method, which assumes that we can approximate the 
signal shape between adjacent samples of the fADC linearly, a 
valid approximation on the leading edge of the pulse. We de- 
termined fo by interpolating the sample before and the sample 
after the halfway point of the maximum sample. Generally, the 
incoming electron in the beam test illuminated several blocks 
simultaneously, allowing us to compute the difference in Iq for 
several different modules. We define the timing difference be- 
tween two modules as Afo.ij - fo,/ - fo,>, where fo^ is the signal 
arrival time Iq for the jth module. 

In principle we can compute the timing difference between 
any two modules, but since we are interested in determining the 
resolution for two modules with comparable signal amplitudes, 
we restricted ourselves to measurements of differences between 
adjacent modules, modules that are either side by side or diag- 
onally aligned with the corners touching. Figure 9 shows the 
distribution of the timing difference Af()_,j for a single module 
and four of its adjacent modules. Previous studies showed that 
the timing resolution depends strongly on the maximum pulse 
amplitude S p, so the data were binned by requiring that for a 
given event, S p for both modules were within the same range. 
For example in Fig. 9, the range of S p was restricted to be be- 
tween 1000 and 2000 ADC counts for the two modules. 
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Figure 9: Distribution of Afo,,) for a single module and four of 
its adjacent modules when all modules had 1000 < S p < 2000 
ADC counts, together with Gaussian fit curves. 



4.2. Determination of timing resolution 

For each range of S p, the distribution of Afo^ was deter- 
mined for all adjacent module combinations. Fits to each distri- 
bution were done with Gaussian functions, such as those shown 
in Fig. 9, to characterize each Afo,,/, which yield the resolu- 
tion in timing difference ctjj with statistical error 5crij. Fig. 10 
shows how the ctij between a given module and four of its ad- 
jacent modules changes with S p. Clearly the timing resolution 
improves with larger pulse amplitude. 
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Figure 10: a-jj between an example module and four of its adja- 
cent modules, determined by fits such as those shown in Fig. 9. 
The four different colors and symbols represent different com- 
binations of modules. Note the log scale on the .^c-axis. 

The value 0",^ is the resolution on the difference in absolute 
time between two modules and is therefore worse than the time 
resolution for a single module. We assumed that the resolution 
of each single module ctj contributes in quadrature to the ob- 
served timing difference resolution so that cr^^. - o-j + cr^., and 
determined each module's resolution from all of the different 
combinations measured. To do this, we divided our modules 
into 2x2 clusters where all modules within the cluster were 
adjacent to each other Within a single cluster there are six pos- 
sible separate measurements of timing differences. If all six 
combinations each had sufficient statistics to yield a fit result. 



we did a fit to minimize 
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where the fit variables cr, gave the intrinsic timing resolution 
of module /. This allowed a determination of the resolution of 
each module within that cluster 

Figure 1 1 shows the resolution of a single module. Because 
each module can be contained in up to four separate clusters, 
the timing resolution in a given bin of S p can be determined up 
to four times. When there were multiple measurements coming 
from different clusters at a given range of S p, the different mea- 
surements were merged by a weighted average, with the stan- 
dard deviation taken as the error. We parametrized the results 
with fits of the form 
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+ b. 



(3) 



where a (ns ■ ADC counts) and b (ns) are fit parameters. Fits 
were also tried with other forms, such as a/ - 



jS p + b and 



a/S p® b, where in the latter form the two terms were added in 
quadrature, but neither of these fit the data better than Eq. (3). 



o 







signal size (ADC counts) 



Figure 11: Final timing resolution for one module. The solid 
magenta points are the weighted averages and standard devia- 
tions of up to four separate clusters that determined the resolu- 
tion for this module. The solid curve shows the fit result, while 
the dotted curve is from Ref. [13]. 

By fitting the timing resolution as a function of S p, we deter- 
mined the parameters a and b for each module. Table 2 shows 
the results of the parameters a and b for the six modules we 
used. By taking the average and standard deviation, our final re- 
sults were a = 72.9±5.3 (ns- ADC counts), b = 0.29+0.02 (ns). 
We can convert the ADC counts into pulse heights in mV: the 
fADCs used 12 bits to digitize the 0.5 V fADC full scale. With 
this conversion, the fit parameter a becomes a = 8.90 ± 0.65 
ns ■ mV. We note that, on average, the gain of the modules is 
such that a signal pulse amplitude of 3 ADC counts corresponds 
to 1 MeV of deposited energy in a block. 

The results from Ref. [13] gave the fit parameters a - 114 ± 
46 (ns ■ ADC counts), b = 0.155 ± 0.077 (ns). The statisti- 
cal parameter can be converted to pulse height with the 1 .45 V 
fADC full scale used, and gives a = 40.4 ± 16.3 ns ■ mV. For 



pulse heights of 100 mV and 500 mV, the timing resolutions 
are then 0.57 + 0.18 ns, and 0.24 + 0.08 ns, respectively. In 
the current measurements the actual pulse heights were limited 
to less than approximately 250 mV, but extrapolating our fit re- 
sults gives timing resolutions of 0.38+0.03 ns and 0.30+0.02 ns 
for 100 mV and 500 mV pulses, respectively. Our results show 
that the timing resolution exceeded expectations at smaller sig- 
nal amplitudes where previous results do not exist. 

Table 2: Measurements of timing resolution parameters a and 
b from Eq. (3). The final row shows the average values and 
standard deviations of the six modules used. 



Module number 


a (ns ■ ADC counts) 


bins) 


7 


77.9 + 8.5 


0.27 ± 0.02 


8 


65.6 + 6.7 


0.32 ± 0.02 


9 


78.9 ±8.1 


0.31 ±0.02 


12 


68.1 ±6.8 


0.31 ±0.02 


13 


72.5 ±7.7 


0.28 ± 0.02 


14 


74.6 ± 7.8 


0.28 ± 0.02 


Average 


72.9 ± 5.3 


0.29 ± 0.02 



In the GlueX experiment, a single photon can cause show- 
ers in multiple blocks, which allow independent timing mea- 
surements from each block. Combining these measurements 
will further enhance the timing resolution of photon showers, 
thereby enhancing our ability to choose which beam bunch pro- 
duced the signal photon of interest and reject background elec- 
tromagnetic showers. 

5. Conclusions 

Results of a beam test using a small version of the FCAL of 
the GlueX experiment are given. The energy resolution is be- 
tween 14% and 20% at energies of 260 MeV down to 1 10 MeV. 
This represents the first test of the production calorimeter hard- 
ware and electronics with actual electromagnetic showers. The 
energy resolution is consistent with the design goals for GlueX 
at these low energies. Precision timing measurements of the 
FCAL signals is a feature that is possible due to the 250 MHz 
fADCs to be used in GlueX. The timing resolution is found to 
exceed the expectations of previous measurements at low PMT 
pulse amplitudes, with timing resolutions of 0.38 ns and bet- 
ter achievable from a single module with signals larger than 
100 mV. A weighted average of timing measurements from dif- 
ferent modules that are illuminated by a single shower will al- 
low enhanced precision of shower time allowing a determina- 
tion of the beam bunch and a rejection of out-of-time electro- 
magnetic background. 
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